Counting atoms using interaction blockade 
in an optical superlattice 
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We report on the observation of an interaction blockade effect for ultracold atoms in optical 
lattices, analogous to Coulomb blockade observed in mesoscopic solid state systems. When the 
lattice sites are converted into biased double wells, we detect a discrete set of steps in the well 
population for increasing bias potentials. These correspond to tunneling resonances where the atom 
number on each side of the barrier changes one by one. This allows us to count and control the 
number of atoms within a given well. By evaluating the amplitude of the different plateaus, we can 
fully determine the number distribution of the atoms in the lattice, which we demonstrate for the 
case of a superfluid and Mott insulating regime of ^^Rb. 
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Interactions between individual particles can drasti- 
cally alter the ground state properties of a many-body 
system and are known to crucially determine the trans- 
port properties in mesoscopic solid state systems. For 
instance, in quantum dots, Coulomb blockade [H El [3l |4] 
can prevent current flow at small bias voltages and in- 
duces tunneling resonances involving a single electron 
crossing the dot. It demonstrates the quantization of 
charge in a striking manner and enables counting and 
control of the charge carriers in the dot with single elec- 
tron resolution, provided that the thermal energy is lower 
than the interaction energy. For the case of ultracold 
atoms, similar interaction blockade effects have been pre- 
dicted [SI [H [71 [HI |9] that could furthermore enable the 
generation of non-classical matter wave states or the 
preparation of highly correlated and entangled samples of 
atoms. Here we present a direct measurement of such an 
interaction blockade with atoms in a bichromatic optical 
superlattice [lOl jTl] [^ . By adiabatically changing each 
site of a tetragonal lattice into a biased double well, the 
atoms experience an interaction blockade forcing them to 
distribute over the two sides of the double wells in a de- 
fined manner. Such an "interacting beam splitter" (IBS) 
for atoms can be used to obtain "constructed pairs" of 
atoms [13^ when slowly recombining a double well, or 
to avoid creating unwanted Fock states when splitting 
[T4|. Observing the interaction blockade-induced tun- 
neling resonances with an intra double well site sensitive 
detection [121 US] provides a novel method of measuring 
the number distribution in the lattice [151 [El [II] • We 
demonstrate this by fully determining the atom number 
distribution of an ultracold degenerate cloud of ^ ''Rb in 
the superfluid or strongly interacting regime of a Mott 
insulator [18 . 

We consider the transport of ultracold atoms between 
the two sides of a double well generated with an optical 
superlattice combined with two transverse lattices [12j . 



The superlattice potential is given by: 

V{x) = Vs cos^(47rx/A^ - 0) + cos^(27rx/A0 , (1) 

where (j) is the relative phase between the short and 
long period lattices. Vs^i^t denote the lattice depths of the 
short, long and transverse lattices, expressed in units of 
their corresponding recoil energies Ep^^^ = /i^/2MA^^^. 
The optical wavelengths are Xs,i,t and M is the atomic 
mass. We consider isolated double wells, which we de- 
scribe with a two-mode Bose-Hubbard type Hamiltonian: 

H = - J (a^^an + aj^a^^ - ^(^l - tir) 

+ ^ {fiL {riL -l)^nR {fiR - 1)) , (2) 

where the operators a]^ ^ and cll^r create and annihi- 
late an atom in the left and right well, ul^r count the 
atom number per well, J is the tunneling matrix element, 
A the potential bias along the double well axis and U the 
onsite interaction energy between two atoms. The state 
of the system can be described in the basis of the Fock 
states Iul^tir), where til^r represent the discrete atom 
numbers on the left and right hand side of the double 
well with a total atom number n = + ur. 

Let us focus on the blockade and calculate for the 
ground state with n atoms |^^)n the probability of each 
atom to be in the left well: {hL)n/^ = n(^^|^L|^^)n/^- 
For non-interacting particles, all Fock states are degen- 
erate at A = and all atoms would gather in the lower 
energy well for |A| ^ J. For U 7^ 0, we diagonalize 
the Hamiltonian ([2| to obtain the eigenstates and corre- 
sponding eigenenergies of the system. For instance, the 
eigenenergies for n = 4 are plotted against the bias A 
in Fig. ^ for U/J = 30. At small bias |A| < /7, the 
ground state of the system mainly consists of |2,2), and 
the transport of an additional atom to the lower well is 
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FIG. 1: (a) Eigenenergies in a double well filled with 4 
atoms versus A for [//J = 30. The single particle resonances 
(large circles) take place between the two energetically lowest 
states. Higher order tunneling resonances occur with smaller 
couplings (smaller circles). Away from these resonances, the 
eigenstates are Fock states, (b) Left-well normalized popula- 
tion versus A for the ground state with 4 atoms (solid line) 
and 3 atoms (dashed line). 



blocked. When A = [/, a tunneling resonance occurs be- 
tween the two states |2, 2) and |3, 1), leaving the latter as 
the ground state for larger bias until the next resonance 
is reached for A = 3U. 

In the experiment, several thousands of double wells 
are populated simultaneously and our measurements can 
only reveal the overall number of atoms in the left (right) 
wells Nl {Nr) and the total atom number N = Nl-\-Nr 
in the lattice. The ratio N^/N is given by the statisti- 
cal average: Nl/N = gn{^L)n/^^ where gn denotes 
the fraction of atoms located in sites with occupation 
n. For a homogeneous system with an integer filling n, 
Nl/N equals the quantum mechanical expectation value 
for a single double well. Such a case is displayed in 
Fig. for n = 3 and n = 4. When the magnitude 
of the bias potential is increased, one observes a succes- 
sion of n steps with an amplitude of 1/n, each corre- 
sponding to one atom changing sides within the double 
well. With i = 1 corresponding to A = 0, we define 

= {Nl[A = {i- 1/2)U] -NL[A = {i- 3/2)U]) /N to 
be the i-th step amplitude, occurring at A = {i — 1)U. 
Even (odd) fillings exhibit only even (odd) steps and the 
step heights for given gn are Si = Y^jL^ 9i+2j / {i + 2j). 
More importantly, the distribution function gn can be 
fully derived from the measured step heights Si as: 



gn = n{Sn- . 



(3) 
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Our experimental sequence begins by loading a ^^Rb 
Bose-Einstein condensate (BEG) into a tetragonal optical 



FIG. 2: (a) Experimental sequence. After loading the BEG 
into a 3D lattice, each site is split slowly into a biased double 
well. The atoms in the left well are transferred into a high 
vibrational level to be counted in a separate Brillouin zone, 
(b) The atoms distribute over the two wells depending on 
the bias A and the interaction energy U. The single particle 
vibrational states are denoted by black lines, whereas states 
shifted by interaction energies 1/2 [/, U etc. per particle are 
shown in red (energy separations not to scale). 



lattice, with periodicities of A//2 = 765 nm on one axis 
and At/2 = 420 nm on the transverse axes. During this 
loading period the lattices are first ramped to = 2 
and Vt = 6 in 100 ms, with all ramps optimized for 
an adiabatic loading [19]. A second ramp over 200 ms 
to final depths of Vi = 40 E^ and Vt = 50 E*^ brings the 
system into a Mott insulating regime. All three tunneling 
frequencies are kept equal during the second ramp, which 
ensures a redistribution of the density on all axes on the 
same time scale. 

The interaction blockade is observed after subsequently 
ramping up the short lattice to Vg = 44 E^ in 20 ms for 
different superlattice phases ^, i.e. different potential bi- 
ases A, as depicted in Fig. [2^. The rising barrier realizes 
an IBS, which forces the atoms to redistribute via tun- 
neling (Fig. ^p) over the two newly-formed wells. 

In order to read out the number of atoms in the left 
well after the IBS sequence, we carry out the detection 
sequence shown in Fig. |2^ (right panel). By increasing 
the bias in 20 ms to a fixed phase of (j)o = — l.lrad, a 
highly tilted double well is formed. Turning off the short 
lattice exponentially with a time constant of 500 /is subse- 
quently transfers the lifted atoms into a high vibrational 
level of the remaining single well [13]. Finally, a band 
mapping technique [T2j [20] allows us to count the popu- 
lations in different energy bands by counting the popu- 
lation of the different Brillouin zones after time-of- flight, 
which directly gives Nl and Nr. The resulting occupa- 
tions Nl/N as a function of the superlattice phase (j) are 
displayed in Fig. [s^i for A/" 2 x 10^. We observe pro- 
nounced plateaus separated by tunneling resonance steps 
up to Ss. For very large bias, one would expect all atoms 
to occupy the lower energy well, however in our measure- 
ments about 5-10% of the atoms are still detected in the 
Brillouin zone corresponding to the side of the double 
well with higher potential energy. We suspect this to be 
a result of excitations created during the fast transfer in 
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FIG. 3: (a) Left well population, measured as a function 
of the phase between the short and long lattices (points), 
theory assuming full adiabaticity during the IBS (dashed line) 
and theory accounting for the dynamical evolution during the 
splitting (solid line). For both theory curves, the amplitudes 
of the steps were extracted from the experimental data. Inset: 
tunneling resonance to higher vibrational levels, (b) Tunnel- 
ing resonances, shown in the derivative of the left well popu- 
lation. Inset: Fock state coefficients in the SPE basis during 
the time evolution (solid) with = 0.06 rad, compared to the 
eigenstates at the corresponding time (dashed). 

the detection step, heating due to spontaneous photon 
scattering in the lattice or collisions during the time-of- 
flight period which can mix the measured populations 
between the Brillouin zones. 

In addition to the resonances induced by the inter- 
action blockade, tunneling resonances to higher excited 
vibrational levels may be observed for very large applied 
potential biases when the system is not in its ground 
state ^21,,22j. In order to reveal their presence, we opti- 
mize the loading to mostly obtain double wells with two 
atoms in the |1,1) Fock state. The left atom is then 
transferred into the first vibrational band of the long lat- 
tice. The IBS and detection sequence are applied again, 
and we obtain the curve displayed in the inset of Fig. [3^ 
which now shows a broad tunneling resonance for a bias 
of A ^ hcoho as expected from theory, when adding the 
second vibrational level to the Fock state basis (solid). 

For the lowest band, we plot the steps expected for 
our lattice configuration and assuming full adiabatic- 
ity during the IBS in Fig. 3a (dashed line). The theo- 
retical result has been obtained from the single-particle 
eigenstates (SPE) of the biased double wells, and using 
as = 5.31 nm for the scattering length [23'. The experi- 
mental and theoretical results clearly show a discrepancy 
in the position and width of the steps, which we attribute 
to a loss of adiabaticity during the splitting process at 
moderately deep lattices. To confirm this, we have per- 



formed a theoretical analysis of the time evolution of the 
splitting process using the basis set of the two lowest 
SPEs. We denote the population of these two states by 
|n^,n — n^)spE7 where Ug is the population of the lower 
energy state. The general state of the n particles can then 
be written as |^(t)) = Ylng=o^ng,n-ng{t)\ng,n-ng)spE' 
We determine the time evolution of the relevant ampli- 
tudes by a numerical solution of the Schrodinger equation 
for the ramps used in the experiment. The resulting the- 
oretical curve is displayed in Fig. [s] (solid line) and the 
inset of Fig. 3b shows the time evolution of the coeffi- 
cients C2,o (black) and Ci^i (red) for a superlattice phase 
of ^ = 0.06 rad close to the S2 resonance. We see that the 
system is not able to completely follow the ground state 
during the splitting process. Both position and width of 
the resonances observed in the experiment are in good 
agreement with the calculation (see Fig. [3] (solid line)) 
when the interaction strength is reduced by 20% from 
the one obtained for the SPE's [24 . We note, however, 
that although the positions of the resonance steps are 
shifted to smaller biases due to this non-adiabatic evolu- 
tion, the amplitude of the steps is fully preserved. The 
relevant atom number distribution is simply frozen out 
at an intermediate short-lattice depth of ^ 20 Ef. during 
the ramp up. 

As outlined in the introduction, the resonance ampli- 
tudes Si can be used to determine the number distribu- 
tion Qn in the initial tetragonal lattice. We have carried 
out such an analysis both for the case of a Mott insu- 
lator and a superfluid as a function of the total atom 
number. We take up to 3 measurements in the center of 
each plateau to obtain the resonance amplitudes Si with 
a typical standard deviation of < 0.5%. In the analysis, 
we have truncated the distribution above a maximum fill- 
ing of n = 4 and rescaled the measured amplitudes by an 
overall factor, such that Si + ^^=2 ~ 

The measured number distribution is displayed in 
Fig. [4^ in the Mott insulator regime, for which we use 
the previous loading sequence. Within our uncertain- 
ties, we detect only an n = 1 component for small 
which corresponds to a single Mott insulator shell. For 
increasing total atom numbers, the next shells appear at 
TV 3 X 10^ for n = 2 and AT 1 X 10^ for n = 3. 
We compare our results with a T = and J = cal- 
culation of the number distribution for our experimental 
parameters |25j and a lattice depth of Vg = ISE^ and 
Vt = 24£^^, which is the configuration where we expect 
the density distribution to be frozen during the ramp up. 
The shaded areas in Fig. [4] correspond to an uncertainty 
of ±10% in the atom number calibration and ±2% in the 
lattice depths and waists. 

In order to measure the number distribution in the 
superfluid regime, we ramp up the long and transverse 
lattices to Vi = GE^ and Vt = 6E^ in a first ramp 
over 200 ms. After a hold time of 40 ms, we freeze 
out the number distribution by ramping the lattices to 
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FIG. 4: Number distribution versus the total atom number, 
measured in the Mott insulator regime (a) and in the super- 
fluid regime (b). The shaded lines are predictions for the 
experimental parameters (see text) and T = 0. Their widths 
represent the uncertainties in atom number and lattice depth 
calibrations. The error bars reflect the standard deviation 
of the measurements. The inset shows a typical occupation 
measurement in the superfluid regime for iV 2 X 10^ with 
four visible steps. 



Vi = 40 El and Vt = 50 in 1 ms. The inset in Fig. ^ 
shows a left- well occupation measurement realized in the 
superfluid regime with ^ 2 x 10^ atoms. Up to the 
fourth resonance step can now be identified. We observe 
that even for the smallest atom numbers, the superfluid 
contains visible fractions of up to n = 3 fillings. A sim- 
ple model which assumes both T = and a Thomas- 
Fermi profile with a Poissonian onsite number distribu- 
tion overestimates the n = 3 and n = 4 fractions. There- 
fore we consider the Gutzwiller approximation adapted 
for bosons [26]. It uses a modified Poissonian distri- 
bution 5^^^~^^A^/n!, which accounts for a possible in- 
teraction induced number squeezing through a squeez- 
ing parameter s = 1/(1 + where z is the number 
of nearest neighbors. For anisotropic lattices, the term 
zJ is replaced by the sum over all nearest-neighbor cou- 
plings. We find that the Gutzwiller distribution models 
our measurement quantitatively for a squeezing parame- 
ter s = 0.66(5) (see Fig.[4]3 shaded areas). The squeezing 
beyond the expected value of s = 0.84 could be due to 
the finite ramp time during the freeze out, or due to 



stronger correlation effects that are not captured by the 
Gutzwiller approximation. 

In conclusion, we have shown how interaction blockade 
can be used to count and control the population of ultra- 
cold atoms on a single side of a double well, both with 
single atom resolution. Furthermore, we have demon- 
strated the ability to measure the full number statistics 
in a 3D optical lattice for different quantum phases using 
this blockade effect. By using our method, it is also pos- 
sible to check predicted phase diagrams in superlattices, 
in particular to verify the existence of fractional filling 
loopholes [27]. Furthermore, by exploiting avoided cross- 
ings in higher order tunneling terms, one could create 
several non-classical and entangled number states within 
the double well in a massively parallel way. 
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